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We address parametric-downconversion seeded by multimode pseudo- thermal fields. We show that this pro- 
cess may be used to generate multimode pairwise correlated states with entanglement properties that can be 
tuned by controlling the seed intensities. Multimode pseudo-thermal fields seeded parametric-downconversion 
represents a novel source of correlated states, which allows one to explore the classical-quantum transition in 
pairwise correlations and to realize ghost imaging and ghost diffraction in regimes not yet explored by experi- 
ments. 
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I. INTRODUCTION 

Ghost imaging [1] and ghost diffraction [2] consist in the 
retrieval of an object transmittance pattern or its Fourier trans- 
form, respectively by evaluating a fourth-order correlation 
function at the detection planes between the field that never 
interacted with the object and a correlated one transmitted by 
the object. A general ghost-imaging/diffraction scheme in- 
volves a source of correlated bipartite fields and two propa- 
gation arms usually called test (T) and reference (R). In the 
T-arm, where the object is placed, a bucket (or a point-like) 
detector measures the total light transmitted by it. The R-arm 
contains an optical setup suitable for reconstructing the image 
of the object or its Fourier transform and a position-sensitive 
detector Jl. 

The correlations needed for ghost imaging and ghost 
diffraction may be either quantum, as those shown by entan- 
gled states produced by spontaneous parametric downconver- 
sion (PDC) H or classical, as those present in the fields at 
the output of a beam-splitter fed with a pair of multi-mode 
pseudo-thermal beams fl d, @]. In recent years several au- 
thors discussed analogies and differences between the two 
cases in terms of the achievable visibility and of the optical 
configurations needed for image reconstruction. An history of 
this debate from different point of views may be found in Refs. 
Jll] and references therein. Recently, it has been suggested 

f~l I - ! rTT 

that the entangled nature of the light source O, [8|, |9fl may 
be necessary to satisfy the "back-propagating" thin-lens equa- 
tion, which, indeed, is fulfilled by PDC-based ghost imaging 
systems. Among other things, we prove that this claim is in- 
correct. 

In this paper, we discuss the use of a novel, PDC-based, 
light source for ghost-imaging/diffraction. In our scheme (see 
Fig. [TJ, the nonlinear crystal realizing PDC is seeded by two 
multi-mode thermal (MMT) beams. We show that the entan- 
glement properties and the amount of correlation at the output 
may be tuned by changing the intensities of the seeds, thus 
leading to a source that can be used to investigate the tran- 
sition from the classical to the quantum regime. Besides, our 



novel source allows ghost-image reconstruction with the same 
optical scheme used for ghost imaging based on spontaneous 
PDC, with the "back-propagating" thin-lens equation that is 
satisfied irrespective of the entanglement of the state. We no- 
tice that the effectiveness of the setup discussed here has been 
already demonstrated in the case of a crystal seeded with a 
single MMT beamlfloll. 

The paper is structured as follows. In section HT1 we calcu- 
late the state obtained from our PDC source with the injection 
of MMT seeds on both T- and R- arms, thus revealing that the 
output field on each arm maintains the statistics of the seed. 
In section [TTTJ we analyze both the intensity correlations be- 
tween the output beams and the entanglement properties of the 
overall state. We explicitly evaluate separability thresholds in 
terms of the seed intensities, and show that the condition for 
the existence of nonclassical correlations in intensity measure- 
ments subsumes the condition for inseparability, i.e sub-shot- 
noise correlations are a sufficient condition for entanglement 
in our system. We also show that entanglement properties of 
the output field are not affected by losses taking place after the 
PDC interaction. In section [TV] we show that the state gener- 
ated in our scheme satisfies the "back-propagating" thin-lens 
equation independently on the seed intensities, i.e. indepen- 
dently on being entangled or not, and it is suitable for realiz- 
ing ghost-imaging and ghost diffraction experiments. Finally, 
Section[V]closes the paper with some concluding remarks. 



II. PARAMETRIC DOWN-CONVERSION WITH 
THERMAL SEEDS 



The interaction scheme we are going to consider is 
schematically depicted in Fig.Q] It consists in a nonlinear x 
crystal pumped by a monochromatic non-depleted plane-wave 
propagating along the z-axis. The Hamiltonian describing the 
resulting parametric process is given by 

Hi = Jd 2 x J^dzx {2) E p (x,t)a T (yi,t)a R (x,t) + h.c. (1) 
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L being the crystal length and x^ 2 ) the nonlinear suscep- 
tibility. The pump-field may be written as E p (x,t) — 
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FIG. 1 : Schematic diagram of the nonlinear interaction. T and R are 
the Test and Reference arms of the setup. 

We can write the interacting quantum fields as 

aj (x, t) oc a j , clj , v .e i ^ z+ ^- r - ( - n ' +u ^ (j = R, T) 

(2) 

where: Kj lZ = yj K"j — q 2 , being the transverse momen- 
tum, Kj = n.j (Tlj + Vj ) /c, rij the index of refraction, fij the 
selected central frequency in channel j, i/j the frequency dis- 
placement with respect to Qj, and c is the speed of light in the 
vacuum. The commutation relation of the quantum fields are 

[aj,<i,u,a}, , ,] = 5j, r 5^5 uy (j,j' = R,T) (3) 



V.q>'] = 



The evolution of a quantum system induced by the inter- 
action Hamiltonian (HJ is described by the unitary operator 

U = exp(— ih~ x J Hj dt), where 



j Hi dt = i^2 Kq, v aT,q, v aR,-<i-u + h.c, (4) 



where K q <v cx sine [(K p — Kt jZ — Kr jZ )L/2]. To obtain Eq. 
(O we have exploited the conservation of energy at the central 
wavelength H p — Ht+^r obtaining vt = — vr = v, and the 



conservation of transverse momentum q^ = — qn = q. As, 
according to Eq. the extension to the non-monochromatic 
case is, in most of cases, straightforward, in the following 
analysis we will focus on the monochromatic emission at the 
frequencies £1 r and Sly and hence we will drop the subscript 
v from the variables. 

The operator U can be rewritten in terms of the operators 
Sq = {nqCLT.cfiR.-ci+h.c.) as U = exp ^ q Sqj . Accord- 
ing to the commutation relations OJ, we have [5 q , Sqi] = 0, 
and therefore U = (2) q e J,Sq , i.e. the interaction establishes 
pairwise correlations among the modes. 

In our analysis we focus on the case in which both the T- 
and R-arms are seeded with MMT beams 



(5) 



Pi 



/j,qj,q\"l i 



n=0 



where j — R,T and |ti) • q denotes the Fock number basis for 
the mode q of the j-arm. The thermal profile of the input is 
given by 

^■,q(n)=M5: q (i+Mi > qr n - 1 . 

/ij.q being the average photon number per mode. The density 
matrix at the output is given by 



Pout 



l e i5q p T , q ® PR^q e iSq , (6) 



According to [? ], it is possible to "disentangle" e lSq by 
using the two-boson representation of the SU(1,1) algebra as 

e lS « =exp{c q 4, q 4,-q} 

x exp |-?7q(a^ q aT,q + aR ) _ q a Bi _ q + 1 )} 



exp{-C q aT, q ai? ,- q } 



(7) 



where £ q = ~«e tanh(|K q |), ?/ q = ln[cosh|«; q |], and 

e^" = Kq/I«q|- 



Equation (0 implies that 



min{m,n} oc 

e lS<1 |n)T,q ® \m) R -q = ^ ^ C q (m, n, k, I) \n-k + l) T ,q ®\m-k + 1)r- c 

k=0 1=0 



(8) 



with 



C (rank I)- e -^(n+m-2fc+i) y/n\m\(n -k + l)\(m - k + 1)1 . fc 



(9) 



3 



min{m,n} oo 

^ ^ Cq(m, n, fa,h) C q (m,n, fc 2 ,Z 2 )* 

:i,fe 2 =0 h,l 2 =0 



By substituting Eq. ([8]) in Eq. ©, we obtain the output state 

min{r, t , It/ u 

Pout = E p ^ n ) p fl,-q( m ) E 1 

\n - fa + h)T,qT,q(n - k 2 + l 2 \ ® \m - fa + h) Ft.-qR,-q(m - k 2 + l 2 \ 
As expected, the first moments of the photon distribution for each mode are those of a thermal statistics 

(«T, q ) = Tr(p OM ia^ q a Tiq ) = /i T ,q + «-PDC,q(l + MT.q + Mfl.-q) 

(™R,- q ) = Tr(p out a^ _ q a fli _ q ) = ^/j,- q + n PDC , q (l + ^T, q + Mfl,-q) 

((An T , q ) 2 } - (n T , q )((n T , q ) + l) 

({AnR^f) = (nji,_ q }((nji,_ q }+l) 

where (O) = Tr[Op oui ] = Tr[U^OUpi n }, AO = O — (O) and ripoc q = sinh 2 |K q | is the average number of photons due to 
spontaneous PDC. 



(10) 



(11) 



r 



Notice that the case of vacuum inputs, pi n = |0)(0|t <8> 
|0)(0|/f, corresponds to spontaneous downconversion, i.e. to 
the generation of twin-beam, whereas the case of a single 
MMT on one arm and the vacuum on the other, pi n = 
® q ( 1 0) (0 T, q 55 Pn.~q) corresponds to the state considered 
inRef. UJ.' 

In the Heisenberg description, the modes after the interac- 
tion with the crystal are given by 6j. q = U^aj ^U, i.e 

6,, q = ZV,,q + e^V q 4__ q (j,f = R, T, j ± f) (12) 
where Uq — cosh|K q | and V q = sinh|K q | (and obviously 

Uq = U-q, V q = V_ q , and Ifq = ip-q). 



III. ENTANGLEMENT AND INTENSITY CORRELATIONS 

In this section we address intensity correlations and entan- 
glement properties of the beams generated in our scheme. As 
we will see, the amount of nonclassical correlations and en- 
tanglement may be tuned upon changing the intensity of the 
thermal seeds and there exist thresholds for the appearance of 
those nonclassical features. On the other hand, the index of 
total correlations (either classical or quantum) is a monoton- 
ically increasing function of the both the seed and the PDC 
energy. 



A. Entanglement and separability 

The downconversion process is known to provide pairwise 
entanglement between signal and idler beams. In our nota- 
tions the (possibly) entangled mode are aT, q and ai?„_ q . In 
the spontaneous process the output state is entangled for any 
value of the parametric gain (i.e for any value of the crystal 
susceptibility, length . . . ) whereas in the case of a thermally 
seeded crystal the degree of entanglement crucially depends 
on the intensity of the seeds. 



Since thermal states are Gaussian and the PDC Hamilto- 
nian is bilinear in the field modes, the overall output state is 
also Gaussian. Therefore the entanglement properties may be 
evaluated by checking the positivity of the partial transpose 
(PPT condition), which represents a sufficient and necessary 
condition for separability for Gaussian pairwise mode entan- 
glement 111 311 . Gaussian states are completely characterized 
by their covariance matrix. In this context let us introduce the 
"position"(-like) operators X and "momentum"(-like) opera- 
tors Y 



X, 



+ °5,q 



V2 



a i a ~ a \ n 



iV2 

Introducing the vector operator 



(13) 



(14) 



with m = 1,2, ...,oo, from the commutation relations (0) 
gives 

Ka,€/9]=*n a , / 3 (15) 

where fl = (J) m a; © a; and ui is the symplectic matrix 



1 

-1 



(16) 



The covariance matrix V is calculated as V a .p = 
2~ 1 ({A£ cn A£ a }), where {Oi,0 2 } denotes the anti- 
commutator. Uncertainty relation among the position and mo- 
mentum operators impose a constraint on the covariance ma- 
trix, V + ifi > 0, corresponding to the positivity of the state. 
The input-output relations for position and momentum opera- 
tors are calculated according to Eq.s (fT2l . obtaining 

C/tX J - q C/=W q X,- q + V q X,v > _ q (17) 

V = W«,Vmi - V„V,. q (j, j' = R, T, J ± f) . 



'q J 3,<i 



q 1 r 
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Without any loss of generality, in the derivation of Eq.s (fTTT i 
we set yj q = which, in turn, corresponds to a proper choice 
of the phase, or, equivalently to a proper redefinition of the 
operators a Jiq corresponding to a rotation of the phase space. 
From Eq.s dTvT > we calculate the covariance matrix 



V 



0Vc 














v q2 











\ 



(18) 



inequality d23T i introduces a threshold. For instance, if we con- 
sider a MMT seed with the same mean number of photon per 
mode, fi, only when the inequality p? > npDc,q(l + 2/i) is 
satisfied, p OM t is separable. It is noteworthy to observe that if 
the PPT is applied to any other subspaces the inequality ob- 
tained are analogous to Eq. d23l l, and thus the result is the 
same. 



B. Separability and losses 



with 








Cq 








./4q 












Bq 













Bq 



(19) 



where 



Aq = [^(2/i T) , 



UqVq(VT,q " 



1) + V 2 (2 MJi ,_ q 

+ 1)- ' 



V^(2/iT, q 

q + 1)- 



l)]/2, 
1)1/2, 



(20) 



V satisfies the uncertainty relations ensuring the positivity of 

Pout ■ 

In order to check whether and when the state p out is en- 
tangled we apply the PPT criteria for Gaussian entanglement 
IU3I1 . For instance, we apply the positive map £_R.- q ' to the 
state p ou t- £-R,-q>(pout) is the transposition (complex con- 
jugation) only of the subspace 7i_R,- q ' corresponding to the 
mode R, — q'. Simon showed that this corresponds to cal- 
culate the covariance matrix V, where all the matrix-blocks 
V q remain the same excepts the matrix V q / — > V q <. V q < 
is calculated with a sign change in the R, — q' momentum 
variable (YR,_ q ' — * — Yr - q ')> while the other momentum 
and position variables remain unchanged (Xx,q' — > -Xr,q', 
>r, q ' — > ir,q'> an d -^H,-q' - * -^fl,-q')- Thus we obtain 



V q ,= 



•Aq' 

























Bq, 








c q < 








(21) 



where _4 q ', Bq' and C q ' are defined in Eq.s (1201 ). According 
to PPT criteria, the separability of p ou u is guaranteed by the 
positivity of L R - ql (p ou t), i.e. 



-n > o 

2 ~ 



(22) 



Inequality ( 122b corresponds to 

MT,q'Mfl,-q' — «PDC,q'(l 4 



A*r, q ' + M-R-q') > °- (23) 



We observe that the spontaneous PDC corresponds to the sit- 
uation with /iT,q' = Mft.-q' = ? thus p onf is entangled. Also 
the case considered in Ref. 11X011 . a MMT seeded PDC only 
on one arm (i.e., p,R.-q> = 0) is always entangled. On the 
contrary, in the case of MMT seeded PDC on both arms, the 



Her we address the problem of the effect of the losses on 
the separability of the state ( fTQb - In fact the presence of losses, 
e.g. internal reflection or absorption in the nonlinear crystal, 
may modify the quantum properties of the state, in particu- 
lar the transition from entanglement to separability (in the ab- 
sence of losses given by the Eq. d23l)). 

Losses in a quantum channel can be modeled by a beam 
splitter in one port of which the quantum channel is in- 
jected while the vacuum enters the other port. The model 
implies that Gaussian states after interaction are still Gaus- 
sian states due to the bi-linearity of the beam-splitter Hamil- 
tonian. Thus also in the presence of losses, the covariance 
matrix completely describes the quantum state. If we con- 
sider an overall transmission factor r on both channels we 
obtain the covariance matrix V r = tV + (1 — r)l/2. The 
form of the covariance matrix V T is completely analogous 
to Eq. ( fT8l ). where the block matrices V q are substituted 
with the block matrices V T q . V T q has the same structure 
of V q in Eq. ( fT9l ), where Aq, £> q , and C q are substituted by 

A T ,q = {1 + 2 T [«| PT,q + V^PR^q + l)]}/2, £ T , q = 

{1 + 2 r pi fi Ri _q + Vl(p T ,q + l)]}/2, and C r , q = r C q , 
respectively. Thus, following the same line of thought of the 
previous section we obtain the covariance matrix V r , corre- 
sponding to the partial transposition of the state. According 
to PPT separability criteria, the state is separable if and only 



if the inequality V T 
be rewritten as 



ii~2 > is fulfilled. This condition can 



T 2 [p T ,q>PR-q- ~ n P DC,q'(l + MT,q' + M-R-q')] > °- ( 2 4) 

Since Ineq. d24l i is fully equivalent to Ineq. d23l l. we conclude 
that losses do not affect the entanglement properties of the 
state in Eq. ( ITOt . 



C. Intensity correlations 

We now evaluate the pairwise intensity correlations owned 
by the generated beams. In addition, we analyze the con- 
nections between threshold for separability and the threshold 
required to have nonclassical correlations. As we will see 
a state obtained by thermally seeded PDC that exhibits sub 
shot-noise correlations is entangled, whereas the converse is 
not necessarily true. In other words, the existence of nonclas- 
sical intensity correlations is a sufficient condition for entan- 
glement. 
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The normalized index of intensity correlation between a 
pair of modes a,- q and a,-', q ' is defined as 



7j,i'(q,q') = 



id' 



(q,q') 



,/((An Tiq )2)((A^- q ) 2 > 
where the correlation term is given by 



(An 3 - iq An 3 v jq / 



(25) 



(26) 



Upon evaluating the first moments as we did in Eq. (fTTb we 
have, for the pair of modes ar, q and a^,_ q , 

rT,fl(q, -q) = "-PDc,q (1 + Jipocq) (1 + Mr, q + Mfl ,-q) 2 



= a 



(27) 



A nonzero value of T^ r, and hence of Jt,r> indicates the 
presence of correlations between the considered modes. Per- 
fect correlations correspond to jt,r = 1. Note that jt.r is 
an increasing function of npDC> an d does not undergo any 
threshold. In Fig [2] we plot -fT,R (solid lines) as a function 
of npDc m two different conditions, namely [it = and 
A*R 7^ (panel (a)) and = fiR ^ (panel (&)). As ex- 
pected 7T,i? approaches unit irrespectively of the mean values 
of the seeding thermal fields as soon as ripnc becomes rel- 
evant. For large npnc q the index of correlation approaches 
unit as follows 



7T,«(q,-q) 



i _ I ^ T >i 

2 



PR, 



2//T,qM.R,-q 1 



(1 + A*T,q + MH,-q) 



'PDC,q 

(28) 



In the two cases jUT, q = A* S> 1 and /i_R,- q = (or viceversa) 
and /iT.q = M-R -q = M ^ 1 we have, respectively 



7T,fl(q, -q) 

7T,fl(q,-q) 



l - 



l 



l 



l 



(1 + n PDC ,q)npDC,q 2/X 
1 „ / 1 



(29) 



(1 + 2n P DC, q ) i 



O 



I*' 



The nonclassical nature of this pairwise correlation may be 
assessed by the quantity 11411 



NRF t ,r{<\) 



((An T , q ) 2 > + ((An fl) _ q ) 2 ) - 2r T , fl (q, -q) 



(»r,q) + (riR-q) 



(30) 



which is usually referred to as "the noise reduction factor". 
A noise reduction, NRFt,r(<i) < 1, indicates the presence 
of nonclassical correlations. The value N RFt.r{ci) = 1 is 
usually called "shot-noise limit" and corresponds to the case 
of a pair of uncorrected coherent signals. By substituting the 
result for our system, we get 



NRF T , R (q) 



MT,q (1 + MT,q) + Mfl,-q (1 + Mfl ,-q) 



MT.q + MR,-q + 2n PDC ,q (1 

We have NRF T , R (q) < 1 if 



MT,q 



1PDC,i 



1 

> - 



a4, 



/4,-q 



2 1 + /i Ti , 



A*R,-q) 

(31) 



(32) 




1E-3 



1E-4 



t 0.01 




PDC 



FIG. 2: The index of total correlations jt,r (solid line) and the noise 
reduction factor NRFt.r (line plus symbol) as a function of npDC 
for the cases of (a): /it = and /j,r ^ and (b) [Xt = ^r 0. 
The values chosen for the parameters are indicated in the figures. 



which subsumes the separability threshold of Eq. (1231 and in- 
dividuates the same region for /iT, q = /i_R,_ q . Therefore, for 
thermally seeded PDC, sub-shot noise correlations imply en- 
tanglement [ 15]. In Fig.|2]we also plot NRFt.r as a function 
of ripDc f° r the same parameters used for Jt,r- As expected, 
the figure shows that NRFt.r, crosses the shot-noise level at 
different values of npnc that depend on the mean values of 
the thermal seeds, thus confirming the intuition that in order 
to achieve sub-shot noise correlations in the presence of two 
thermal seeds we need to have a PDC process strong enough. 



IV. MMT-PDC BASED GHOST IMAGING AND GHOST 
DIFFRACTION 

The bipartite state obtained by the nonlinear process 
described above is suitable for applications to ghost- 
imaging/diffraction protocols. Ghost-imaging and ghost 
diffraction protocols rely on the capability of retrieving an ob- 
ject transmittance pattern and its Fourier transform, respec- 
tively, by the evaluation of a fourth-order correlation function 
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at the detection planes of a light field that has never interacted 
with the object and a correlated one transmitted by the object. 
We consider the schemes depicted in Fig. [3] An object, de- 
scribed by the transmission function t(x T ), is inserted in the 
T-arm on the plane x T and a bucket detector measures the total 
light, It, transmitted by the transparency. The R-arm contains 
an optical setup suitable for reconstructing either the image 
of the object or its Fourier transform and a position-sensitive 
detector that measures the local I R (x R ). The procedure for 
calculating the correlation function between the light detected 
in the two arms of the setup is equivalent to evaluating first 
the correlation function between Ir(xr) and It(xt): 



G^(xr,x t ) = (A/ fl ( Xil )AJ T (x r )), 
and then integrating over all the values of x-p 

S (2) (x fl ) = J dx T G^(x R ,X T ), 



(33) 



(34) 



where (Jjfo)) = (ctfejcffe)) (j = R,T) is the 
mean intensity of the y'-th beam at the detection plane, with 
hereinafter (...) = Tr(...p in ), and (Ir(xr)I t (x t )) = 
(c R (x R )c R (x fl )c T (x t )ct (x t )) . 

The connection between the field operators at the detection 
planes and those at the output of the crystal is given by 



Cj ( Xi ) = J dx^ hj (xj , x'j ) bj (xj ) , 



(35) 



where bj (xfj ) are the Reference and Test field operators at the 
output face of the crystal and Iir(xr,x r ) and /ix(xx)X^) 
are the two response functions describing the propagation of 
the field in the two arms of the setup 11611 . 



By using Eqs. ( f35l > and d34l ) we can rewrite G^ 2 ' (xr, xt) as 

G (2) (xr, x t ) = J dx R dx R dx' T dx T h R (x R , x' R ) h* R (xR, x R ) h T ix T , x' T ) h T (x T , x T ) 
x \bUx' R )b R (x' R )b T (x' T )b T (x' T )) - (bUx' R )b R (x' R )){bUx' T )b T (x' T )) 



(36) 



Also in this case the factorization rule for (b R (x R )b r (x' b ) b T (x^ )6y (x r )) is the same as that for spontaneous PDC lfl7ll . and 
for multi-thermal one-arm-seeded PDC [ 10] 



(b R (x' R )b R (x' R )bU^T)bT(x' T )) = (b R (x' R )b R (x> R ))(b T (x'T)bT(x> T )) + (b R (x' R )bUx'T))(bR(x' R )b T (x T )). 



(37) 



The result in Eq. d37l > can be demonstrated by rewriting bj (x) in terms plane waves as bj (x) <x e * q x ^j,q an d men exploiting 
the input-output relation of Eq. (Y2\ . According to Eq. (I?7] |. and in complete analogy with the case of spontaneous PDC| 17], 
also in the case of the MMT-seeded PDC we obtain 



where 



G^(x R ,x T ) = I dx' R J dx' T h R (x R ,x' R )h T (x T ,x' T )(b R (x' R )b T (x' T )) 
(bR(x' R )b T (x T )) oc V e »[q-(x^-xy +Vq ] WqVq (1 + + = ^[q-^-^+fJ Cq 



q 



q 



(38) 



(39) 



and C q is calculated in Eq. ( f20b By using Eq. d39l . Eq. (l38l can be rewritten as 



x R ,x T ) 



E 

q 



h R (x R , -q)h T (x T , q)C q 



(40) 



where h 3 (xj , q) = / dx^ e' q x j hj (Xj , x'j ). 

According to Fig. [3] we consider two different schemes for the collection optics in the Test arm of the setup: 

(a) the detection plane coincides with the plane of the transparency, x T — xy, and hence 

Mx T ,q) oc e" ! ^ 92 e' !qXT t(xT) 
only describes free propagation over a distance d\\ 

(b) a collection lens is located on the plane x/ y beyond the transparency in a Fourier-transform configuration, and hence 



2tt 



x t ) 



(41) 



(42) 
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FIG. 3: Experimental setup for ghost imaging: t(x-r), object transmission function; fr,R- (Left): experimental configuration with detection 
plane coinciding with the object plane. (Right): experimental configuration with detection plane coinciding with the Fourier plane of the 
collecting lens in the Test arm. 



The optical scheme in the Reference arm contains a lens (focal length f R ) located at x; >R and thus the Fourier transform of the 
impulse response functions can be written as 

oc [ <bt lR e-i^ R+ci y* l - R e if ^-^)<*. (43) 



If fn, =h d-3, Eq. (|43l becomes 

ftfl(x fl ,-q) cx e -'^ (*»+ 1733^1775 ) «V i/d 3 -i// H i' x « ; ( 44 ) 



while if //j = c? 3 , Eq. d43l becomes 



■ Ad 2 2 / 27T 

h R (x R , -q) oc e 1 c ? <J —x R + q . (45) 



Ad 



Depending on the chosen geometrical configuration, these schemes realize either a ghost-imaging or a ghost-diffraction experi- 
ment HI]. 



A. Ghost imaging 



To perform a ghost-imaging experiment we choose f R ^ d^. Let us first consider case (a). Substituting Eq. d4H and Eq. ( f44b 
into Eq. (f40b yields the expression 



G( 2 )(x K ,x T ) oc |i(x T )| 5 



t(x T )| 2 |C q | 2 J(x T + 



M 



which, once integrated over the bucket detector, 



G^(x R ) = / dx T G( 2 )(x i? ,x T )c t(--^) |C q | 2 , 



(46) 



(47) 
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gives the image of the object. Note that in passing from Eq. d40t to Eq. we have made the following assumptions: C q is 
almost independent of q and the distances d\, d-2, and d 3 satisfy the so-called "back-propagating thin lens equation", l/(d\ + 
^2) + 1/^3 = 1//-R [19], so that we obtain an imaging system with magnification factor M = d^/{d\ + c?2)- In case (6), that 
is with a collection lens in the Test arm we proceed similarly by substituting Eq. d42b and Eq. d44b into Eq. d40b and making the 
same assumptions as before: 



G< 2 Wx T ) 



oc 



|«-*)| W 



(48) 



Note that in this case (b) the image of the object emerges from correlations without need to perform the integration over the 
bucket detector in order to recover the ghost image. 



B. Ghost diffraction 



To perform a ghost-diffraction experiment we consider the configuration d% — and choose configuration (b) in the Test 
arm. By inserting Eq. (l42l and Eq. (l45l l into Eq. (l40b 



G^x^xt) 



2?r 2tt 
Ad 3 A/t 



|C q | 



By selecting the component xj- = on the Test plane we get 



t{ -xd-*«' 



|c q | 2 , 



(49) 



(50) 



which gives the diffraction pattern of the object. Note that the choice (a) would not give any meaningful result. 



V. CONCLUSIONS AND OUTLOOKS 

This paper was aimed at showing the possibility of perform- 
ing ghost imaging and ghost diffraction with a novel source 
based on PDC seeded with two MMT fields, which gener- 
ates a bipartite correlated state. Peculiar properties of this 
new source may open a new insight into the understanding 
of the ghost imaging/diffraction process. In fact, nowadays 
the sources considered for ghost imaging/diffraction either 
were definitely separable (classically correlated beams ob- 
tained from a MMT source) or entangled (spontaneous PDC). 
On the contrary, here we proved that the separable/entangled 
nature of the light produced by our source can be controlled 
by changing the seed intensities and that the transition from 
quantum to classical regimes does not modify the possibility 
of realizing ghost imaging schemes. 

Furthermore, we also showed that a ghost imaging ex- 
periment performed with our source satisfies the"back- 
propagating" thin-lens equation, as much as with spontaneous 



PDC, even when the state produced becomes separable. This 
is in contrast with the idea, also recently suggested J3>|H0], 
that the "back-propagating" thin-lens equation is connected 
with the entangled nature of the spontaneous PDC. 

According to the consideration of above, we are planning to 
realize a ghost imaging experiment with a MMT seeded PDC 
source in order to show that the same optical configuration 
allows retrieving of the image irrespectively to the entangle 
or separable nature of the light produced by the source. This 
will definitely demonstrate that there is not any connection 
between the entangled properties of the light source and the 
"back-propagating" thin-lens equation. 
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